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Abstract. A pleasant family of graphs defined by Godsil and McKay is shown to have 
easily computed eigenvalues in many cases. 

Let G and H be directed graphs on the respective vertices U and V, and suppose that the vertex 
sets have each been partitioned into disjoint subsets U = Uq U U\ and V = Vo U V\. The partitioned 
tensor product GxH of G and H with respect to this partitioning is defined as follows: 

a) Each vertex of Uo is replaced by a copy of H\Vo, the subgraph of H induced by Vo; 

b) Each vertex of U\ is replaced by a copy of H|Vi; 

c) Each arc of G that runs from Uo to U\ is replaced by a copy of the arcs of H that run from Vb 
to V x \ 

d) Each arc of G that runs from Ux to Uq is replaced by a copy of the arcs of H that run from V\ 
to Vq. 



G 



H 



GxH 




(a) 

Figure 1. Partitioned tensor products, directed and undirected. 
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For example, Figure 1 shows two partitioned tensor products. The example in Figure lb is undi- 
rected; this is the special case of a directed graph where each undirected edge corresponds to a pair 
of arcs in opposite directions. Arcs of G that stay within Uq or U\ do not contribute to GxH, so 
we may assume that no such arcs exist (i.e., that G is bipartite). 

Figure 2 shows what happens if we interchange the roles of Uq and U\ in G but leave everything 
else intact. (Equivalently, we could interchange the roles of Vq and V\.) These graphs, which may 
be denoted G R xH to distinguish them from the graphs GxH in Figure 1, might look quite different 
from their left-right duals, yet it turns out that the characteristic polynomials of GxH and G R xH 
are strongly related. 



Let Eij be the arcs from Ui to Uj in G, and Fij the arcs from Vi to Vj in H; multiple arcs 
are allowed, so and Fij are multisets. It follows that GxH has |J7o| \ Vo\ + \Ui \ |Vi| vertices 
and \U \ \F 00 \ + |*7i| |-Fn| + |^oi| I -Foil + l-^iol |-Fio| arcs. Similarly, G R xH has \Ut\ \V \ + \U \ |V"i| 
vertices and \Ui\ \F 00 \ + \Uq\ \Fu\ + \E 10 \ \F 01 \ + \E 01 \ \F 10 \ arcs. 

The definition of partitioned tensor product is due to Godsil and McKay [3] , who proved the 
remarkable fact that 



where p denotes the characteristic polynomial of a graph. They also observed [4] that Figures lb 
and 2b represent the smallest pair of connected undirected graphs having the same spectrum (the 
same p). The purpose of the present note is to refine their results by showing how to calculate 
p{GxH) explicitly in terms of G and H. 

We can use the symbols G and H to stand for the adjacency matrices as well as for the graphs 
themselves. Thus we have 



G R xH = 




(a) 



Figure 2. Dual products after right-left reflection of G. 



p(GxH)p(H\V )\ Ul \- ]u °\ =p(G T xfl>(#1^i) |C/lHC/o1 , 




and 
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in partitioned form, where G^ and Hij denote the respective adjacency matrices corresponding to 
the arcs Eij and Fij. (These submatrices are not necessarily square; Gij has size \ Ui \ x \Uj \ and Hij 
has size |V^| x \Vj\.) It follows by definition that 

GxH = fVol ® H oo G i®H 01 \ 
~ ~\G 10 ®H 10 I^^Hu) 
where <8> denotes the Kronecker product or tensor product [7, page 8] and Ik denotes an identity 
matrix of size k x k. 

Let H | a denote the graph obtained from H by cr-fold repetition of each arc that joins Vb 
to V\. In matrix form 

This definition applies to the adjacency matrix when a is any complex number, but of course 
H | a is difficult to "draw" unless cr is a nonnegative integer. We will show that the characteristic 
polynomial of GxH factors into characteristic polynomials of graphs H f cr, times a power of the 
characteristic polynomials of H 00 or H\\. The proof is simplest when G is undirected. 

Theorem 1 . Let G be an undirected graph, and let (<7i , . . . , cr; ) be the singular values of Goi = Gf , 
where Z = min(|J7 |,|?7i|). Then 

f(nU^T^))p(^oo) l%| - |c/l1 , if |£/o| > l^il; 

p(Gxtf) = { ) ( 

[(nJ=iP(^T^))p(^ii) |t7i| - |I/o1 , if i^ii>i^oi. 

Proof. Any real mx n matrix A has a singular value decomposition 

A = QSR T 

where Q is an m x m orthogonal matrix, R is an n x n orthogonal matrix, and S is an m x n 
matrix with Sjj = Oj > for 1 < j < min(m, n) and Sij = for i / j [6, page 16]. The numbers 
cti, . . . , cr min ( m n ) are called the singular values of A. 

Let m = \Uo\ and n = and suppose that QSR T is the singular value decomposition of Goi- 
Then . . . , gi) are the nonnegative eigenvalues of the bipartite graph G, and we have 

m T ®I\ Vo \ O \ GxH (Q® 1 ^ ° \ = (I\u \ ®Hqo S®H 01 \ 

V o i? T ®i| Vl |y - v o R®i\ Vl \) \s t ®h w i^qhuJ 

because Gio = RS T Q T . Row and column permutations of this matrix transform it into the block 
diagonal form 

/HUi \ 

\ Dj 

where D consists of m — n copies of Hqq if m > n, or n — m copies of H\\ if n > m. □ 

A similar result holds when G is directed, but we cannot use the singular value decomposition 
because the eigenvalues of G might not be real and the elementary divisors of XI — G might not 
be linear. The following lemma can be used in place of the singular value decomposition in such 
cases. 
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Lemma. Let A and B be arbitrary matrices of complex numbers, where A is m x n and B is 
n x m. Then we can write 

A = QSR- 1 , B = RTQ- 1 , 

where Q is a nonsingular m x m matrix, R is a nonsingular n x n matrix, S is an m x n matrix, 
T is an n x m matrix, and the matrices (S, T) are triangular with consistent diagonals: 

Sij = Tij =0 fori>j; 

Sjj = Tjj or SjjTjj = , for 1 < j < min(m, n) . 

Proof. We may assume that m < n. If AB has a nonzero eigenvalue A, let a be any square 
root of A and let x be a nonzero m- vector such that ABx = a 2 x. Then the n- vector y = Bxjo is 
nonzero, and we have 

Ay = ax , Bx = ay . 

On the other hand, if all eigenvalues of AB are zero, let x be a nonzero vector such that ABx = 0. 
Then if Bx / 0, let y = Bx. If Bx = 0, let y be any nonzero vector such that Ay = 0; this is 
possible unless all n columns of A are linearly independent, in which case we must have m = n 
and we can find y such that Ay = x. In all cases we have therefore demonstrated the existence of 
nonzero vectors x and y such that 

Ay = ax , Bx = ry , a = t or ar = . 

Let X be a nonsingular m x m matrix whose first column is x, and let Y be a nonsingular 
n x n matrix whose first column is y. Then 

where A± is (m — 1) x (n — 1) and B 1 is (n — 1) x (m — 1). If m = 1, let Q = X, R = Y, S = (aa), 
and T = Qj). Otherwise we have A\ = QiSiR^ 1 and Bi = RiTiQ^ 1 by induction, and we can let 

£)• «-(o £)• s Ko £)• r =(; B *)- 

All conditions are now fulfilled. □ 

Theorem 2. Let G be an arbitrary graph, and let (a\, . . . ,07) be such that aj = Sjj = Tjj or 
a j = = SjjTjj when Gqi = QSR -1 andGio = RTQ -1 as in the lemma, wherel = min(|C/o| , \U\\). 
Then p(GxH) satisfies the identities of Theorem 1. 

Proof. Proceeding as in the proof of Theorem 1, we have 

fQ- 1 ®^ O \ GxH (Q® I Wo\ O \ = n\ Uo \®Hoo S®H 01 \ 

V O R- 1 ®!^) ~ V O R®I\ Vl \) V T®H W I\ Vl \ ®Hu) ' 
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This time a row and column permutation converts the right-hand matrix to a block triangular 
form, with zeroes below the diagonal blocks. Each block on the diagonal is either H f Oj or Hqq 
or Hn, or of the form 

( H 00 a H i 



\tHW H n " aT = °- 
In the latter case the characteristic polynomial is clearly p(H 00 )p(Hu) = p{H | 0), so the remainder 
of the proof of Theorem 1 carries over in general. □ 

The proof of the lemma shows that the numbers a 2 ,..., a 2 , are the characteristic roots of 
GoiGiO; when \Uo\ < \Ui\, otherwise they are the characteristic roots of GioGoi- Either square 
root of (Tj can be chosen, since the matrix H | a is similar to H f {—&)■ 

We have now reduced the problem of computing p(G x H) to the problem of computing the 
characteristic polynomial of the graphs H f a. The latter is easy when a = 0, and some graphs G 
have only a few nonzero singular values. For example, if G is the complete bipartite graph having 
parts Uq and U\ of sizes m and n, all singular values vanish except for a = ^Jmn. 

If H is small, and if only a few nonzero a need to be considered, the computation of p(H f cr) 
can be carried out directly. For example, it turns out that 
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A/ 



(A 2 + A - a 2 )(X 3 - A 2 - (2 + a 2 )\ + 2) ; 



so we can compute the spectrum of G x H by solving a few quadratic and cubic equations, when 
H is this particular 5- vertex graph (a partitioned 5-cycle). But it is interesting to look for large 
families of graphs for which simple formulas yield p{H f a) as a function of a. 

One such family consists of graphs that have only one edge crossing the partition. Let H 00 
and Hn be graphs on Vq and V\, and form the graph H = Hqq — H\\ by adding a single edge 
between designated vertices xo € Vq and x\ G V\. Then a glance at the adjacency matrix of H 
shows that 

p(H T a) = p(fToo)p(#n) - aVtfoolVbWMtfiilVAx!) . 

(The special case cr = 1 of this formula is Theorem 4.2(h) of [5].) 

Another case where p(H f a) has a simple form arises when the matrices 

#oo A , „ _ / #oi 



Fo= l o Hll ) and ^"V^io o 

commute with each other. Then it is well known [2] that the eigenvalues of Hq + (tH± are A,- +afij, 
for some ordering of the eigenvalues Xj of H and fj,j of Hi. Let us say that (Vq, Vi) is a compatible 
partition of i7 if H Hi = HiH 0l i.e., if 

HqqHqi = HqiHu and HnH w = H W H 00 . 

When if is undirected, so that #00 = #00 an d #11 = #11 an d #10 = #oi > the compatibility 
condition boils down to the single relation 

HqqHqi = HqiHh . (*) 
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Let m = | Vb| and n = \V\\, so that Hqq is m x m, H 0i ismxn, and Hu is n x n. One obvious 
way to satisfy (*) is to let Hqq and Hu both be zero, so that H is bipartite as well as G. Then 
H | a is simply u-ff , the c-fold repetition of the arcs of i?, and its eigenvalues are just those of H 
multiplied by a. For example, if G is the M-cube and H is a path P^f on iV points, and if Uq 
consists of the vertices of even parity in G while Vq is one of -ff's bipartite parts, the characteristic 
polynomial of G x H is 

n ^-(^-^^ ■ 

l<j<M V 7 
l<fc<JV 

because of the well-known eigenvalues of G and H [1]. Figure 3 illustrates this construction in the 
special case M = N = 3. The smallest pair of cospectral graphs, and p~|, is obtained in a 
similar way by considering the eigenvalues of P3 x P 3 and Pj x P3 [4] . 




Another simple way to satisfy the compatibility condition (*) with symmetric matrices Hqq 
and Hu is to let Hqi consist entirely of Is, and to let Hqq and Hu both be regular graphs of 
the same degree d. Then the eigenvalues of H are (Ai, . . . , A m , X[, . . . , A^), where (Ai, . . . , X m ) 
belong to Hqq and (A2,...,A^) belong to Hu and Ai = A' : = d. The eigenvalues of Hi are 
{yjmn, —yjmn, 0, . . . , 0). We can match the eigenvalues of H properly with those of Hi by looking 
at the common eigenvectors (1, . . . , 1) T and (1, . . . , 1, —1, . . . , — 1) T that correspond to d in H and 
±^/mn in Hi ; the eigenvalues of H f a are therefore 

(d + ay/mn, A2, . . . , A m , d — ay/mn, X' 2 , ■ ■ ■ , A^) . 

Yet another easy way to satisfy (*) is to assume that m = n and to let H 00 = Hu commute 
with Hq\. One general construction of this kind arises when the vertices of Vq and Vi are the 
elements of a group, and when Hoo = Hu is a Cayley graph on that group. In other words, 
two elements a and (3 are adjacent in Hqq iff a(3~ l G X, where X is an arbitrary set of group 
elements closed under inverses. And we can let a £ Vq be adjacent to (3 G Vi iff a/3 -1 G Y", 
where Y" is any normal subgroup. Then ff o commutes with i7 i- The effect is to make the 
cosets of Y fully interconnected between Vq and Vi, while retaining a more interesting Cayley 
graph structure inside Vq and V\. If Y is the trivial subgroup, so that H 01 is simply the identity 
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matrix, our partitioned tensor product G xH becomes simply the ordinary Cartesian product 
G © H = (g) H + G (g) I\y\. But in many other cases this construction gives something more 
general. 

A fourth family of compatible partitions is illustrated by the following graph H in which m = 6 
and n = 12: 
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In general, let C2k be the matrix of a cyclic permutation on 2k elements, and let m = 2k, n = 4k. 
Then we obtain a compatible partition if 



Hqq — {C 3 2k + C\ k + C 2fc J ) , H i — (I 2 k Cik) , Hu — ^ 



r<3 i ri-j /~ik+l 



2fc 2/c 2/c 

2k + v .7, 



(~tk-\ (~<j , (-t-j 



The 18 x 18 example matrix is the special case j = 2, k = 3. The eigenvalues of H f a in general 
are 

u jl + uj~ jl + 1 , uj jl + uj- j1 -1 + V2a, J 1 + uj~ jl -l-V2a 

for < I < 2k, where u = e ni / k . 

Compatible partitionings of digraphs are not difficult to construct. But it would be interesting 
to find further examples of undirected graphs, without multiple edges, that have a compatible 
partition. 
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